Ferromagnetism in the Strong Hybridization Regime of the Periodic Anderson Model 
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We determine exactly the ground state of the one-dimensional periodic Anderson model (PAM) 
in the strong hybridization regime. In this regime, the low energy sector of the PAM maps into an 
effective Hamiltonian that has a ferromagnetic ground state for any electron density between half 
and three quarters filling. This rigorous result proves the existence of a new magnetic state that 
was excluded in the previous analysis of the mixed valence systems. 



Rigorous results, numerical and analytic, have greatly 
aided the study of strongly correlated electrons systems. 
Unfortunately, few such results exist. The numerical 
renormalization group and Bethe ansatz solutions of the 
single impurity Anderson and Kondo models are perhaps 
the best examples of a solid numerical approach and ex- 
act solution of simple models that changed and solidified 
the thinking in what was a highly controversial and puz- 
zling problem area. Another important result was the 
rigorous connection between the two models established 
by Schrieffer and Wolff [1]. Their work showed that the 
low lying energy spectrum of the Anderson model in the 
strong coupling and weak hybridization limit (U/t S> 1 
and V <C \Ep—ef\) can be mapped into the Kondo model 
in the weak coupling regime (J ft <C 1). 

For dense systems, the natural extensions of the impu- 
rity models are the periodic Anderson (PAM) and Kondo 
lattice (KLM) models. For these, numerical renormaliza- 
tion group and Bethe Ansatz solutions are lacking even 
in one-dimension. In addition, very few rigorous results 
are available for the PAM [2,3]. What remains true how- 
ever is the connection between the strong and weak cou- 
pling limits via a natural extension of the Schrieffer- Wolff 
transformation. 

There are two basic questions one can ask about these 
lattice models: what is their relevance to real materials 
and in what other parameter regimes might their physics 
be connected? In several well known papers, Doniach [4], 
at least implicitly, made several assumptions about the 
answers to both questions and proposed the now stan- 
dard picture of the magnetic properties of f-electron ma- 
terials that portrays a competition between the RKKY 
magnetic interaction which is obtained from a fourth or- 
der expansion in the hybridization and the Kondo ex- 
change. 

The reasoning behind this intuitively appealing picture 
is something like the following: From the Schrieffer- Wolff 
perturbation theory, the Kondo exchange coupling is re- 
lated to the parameters in the PAM via Jsa \V\ 2 /\Ep — 
ef\. In the PAM, a mixed valence regime corresponds 
to positioning the f-electron orbitals in the conduction 
band near the Fermi energy, i.e., \Ep — €f \ ~ 0. In this 



regime, the Schrieffer- Wolff result suggests the Kondo 
exchange is strong and thus can lead to a complete com- 
pensation of the f-moments by one or more of the con- 
duction band electrons. Implied in this line of reasoning 
there are two important assumptions [4]: I) The strong 
coupling limit of the KLM is connected to the mixed va- 
lence regime (for weak hybridzation \V\ -C \t\) of the 
PAM, and II) the number of conduction electrons per 
/-magnetic moment is larger than one. Regarding the 
first assumption, we note that this connection is not es- 
tablished by the Schrieffer- Wolff transformation because 
when \Ep — e/| w 0, the transformation is no longer valid. 

We also note recent numerical studies [5,6] of the weak 
hybridzation, strongly coupled PAM that find over a wide 
range of parameters ferromagnetic states which are not a 
result of the RKKY interaction and non-magnetic states 
which are not a consequence of a Kondo-like compen- 
sation of the f-moments by the conduction band elec- 
trons [7]. Additionally, there are mixed valence materi- 
als which exhibit a co-existence of ferromagnetism and 
a strong Kondo-like behavior (e.g., CeSLj [8], CeGe2 [9], 
CetRh^Ru^Ba [10], CeSii. 76 Cu . 24 [11], Ce 3 Bi 4 , and 
CeNio.sPto.2 [12]). Because the RKKY interaction is 
considerably weaker than the Kondo exchange, these ex- 
perimental results suggest that a ferromagnetic mecha- 
nism other than RKKY, at least sometimes, dominates 
in mixed valence materials. 

In this paper we note that a strong hybrization, |V| 3> 
\t\, can lead to mixed valence state and a one-on-one 
compensation of an f-moment by a conduction electron. 
Then, for the strong couping limit of the one-dimensional 
KLM we rigorously establish that the ground state is con- 
nected to the strong hybridzation, strong coupling limit 
of the PAM. We thus show that Doniach's first assump- 
tion is valid if the hybridization is strong, although its va- 
lidity remains questionable for the more relevant regime 
\V\ <C \t\ (weak hiybridization) . In the strong hybridiza- 
tion limit we rigorously show that there is ferromagnetic 
ordering instead of a nonmagnetic Kondo state. The 
difference between our result and Doniach's can be at- 
tributed to the violation of his second assumption which 
indeed is not valid for a lattice system (one /-orbital per 
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unit cell) since the number of conduction electrons (with 
net magnetic moment) per / spin cannot be larger than 
one (it is equal to one only at half-filling) . 

The Hamiltonian for the one dimensional PAM is: 
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have energy = e/. 

For three particles on one site, there are two possible 
states due to the two possible orientations of the / spin: 
fi a d\^d\^\0) . The energy of both states is E 3 — Cf. 



where d\ a and create an electron with spin a in the 
d and f orbitals of the lattice site i and n{„ = f,„ f-. 

•> 1(7 ^ 1(7 " 1(7 

We will show that for infinite U, \t\ <C |V|, and 
l e /l ^ \y\ (asymmetric regime), the low energy sector 
of the PAM can be mapped into an infinite U Hubbard 
model which includes correlated next-nearest neighbor 
hoppings and nearest-neighbor repulsions. To this end 
we first need to solve the atomic (one site) limit of the 
PAM [13] for all the possible fillings (0 to 3 particles per 
site because U is infinite). 

There are two possible eigenstates for one particle on 
one site i. These eigenstates are created by the following 
operators: 
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The operator a[ a creates a particle in the bonding state 



with energy E 1 = ^ — y ^ + V 2 , while (3\ a creates 
a particle in the anti-bonding state with energy E^ = 



T + Vf + V2 ( see Fi §' 

For two particles on one site, there are two singlets 

(U = oo), and three triplet eigenstates (see Fig. 1). The 

ground and the highest energy states are the bonding and 

the anti-bonding singlets respectively: 
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The energies of the bonding and the anti-bonding singlets 



are: E^ 
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+ 2V 2 . The triplet states: 
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FIG. 1. Energy levels for the atomic solution of the PAM 
for infinite U. The full lines correspond to the lowest energy 
levels that generate the reduced subspace of H e ff. 

We will consider the range of concentrations: \ < n < 
|, where n = N e /AL and N e is the total number of elec- 
trons. In this way, the concentration ranges from two to 
three particles per site. For t = 0, the ground state of H 
is massively degenerate, and the corresponding subspace 
is generated by states containing L(3 — An) local bond- 
ing singlets \cf)~) (doubly occupied sites) and 4L(n — i) 
sites occupied by three particles. Therefore, the effec- 
tive Hamiltonian which is obtained when t is included 
perturbatively has a local dimension equal to 3 because 
each site can be occupied by a bonding singlet (hole) 
or by three particles state with two possible spin orien- 
tations (S=l/2 particle). The huge degeneracy of the 
ground state is then associated with both the spin and 
the charge degrees of freedom. We will see below that 
while the degeneracy associated with the charge degrees 
of freedom is lifted to first order in t, the spin degeneracy 
is lifted at second order. 

To derive an effective Hamiltonian, we first need to 
identify the possible virtual processes. There are two 
different types of virtual (excited) states which are ob- 
tained by applying the hopping term to the ground state 
subspace of H(t — 0) (see Fig. 2): a) those in which 
two local nearest neighbor bonding singlets (2 + 2) are 
excited to local states having one and three particles 
(1 + 3) and b) those in which a local bonding singlet and 
a three particle state, which are nearest neighbors, are 
permuted and excited to the other possible local states 
(anti-bonding singlet or triplet). In any of these virtual 
processes, t is much smaller than the energy difference be- 
tween the virtual and the ground state because \t\ <^ \V\ 
and \e f \ < \V\. 
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FIG. 2. Virtual states which are obtained from the appli- 
cation of the hopping term t to the lowest energy subspace of 
H(t = 0). 

Introducing the creation and annihilation operators: 

4 = (i ( 6 ) 



and mapping the bonding singlet, \4>~), into the empty 
state or hole, we can write the effective Hamiltonian up 
to second order in t: 



*1 X^^la^i-lff + #- C 
i,<7 

4 X! ( C l+l<r' c L c i<r' c i-l<r + H - c ) 

i,0^,0^ , 

9^2(1 - «i+i)(l - Tlx) + e f L 

i.a 



(7) 



with 



2^2 



t z a 



2' 1 2(E^-E 2 )' 

2" 



(8) 



^2 ~ -^2 



-^2 ~ E 2 



{au + V2vb) 2 (av - V2ub) 2 



2Eo 



2E; 
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The half-filled PAM corresponds to the absence of c elec- 
trons (vacuum) and n = | corresponds to one c electron 
per site. 



If we only keep the first order terms of H e ff, the re- 
sulting Hamiltonian is an infinite U Hubbard model. In 
this model, the charge and the spin degrees of freedom 
are decoupled and therefore the wave function is a prod- 
uct of a charge and a spin component. As a consequence, 
the infinite U Hubbard Hamiltonian maps into a spinlcss 
model for each fixed spin configuration. Therefore, the 
complete spin degeneracy persists to first order in t while 
the charge degeneracy is lifted. 

The second order terms introduce next nearest neigh- 
bor correlated hoppings t' l7 t' 2 , t' 3 and a repulsion g be- 
tween nearest neighbors. The terms t[ and t' 2 describe 
hopping processes where the electron hops over another 
electron in the intermediate site with and without spin 
flip. The t' 3 term describes the hopping over an empty 
site. Notice that t[ is the only term which can change 
the spin ordering. 

Using a second order expansion around the strong cou- 
pling limit (large J/t) of the Kondo lattice model (KLM), 
Sigrist et al [14] found an effective Hamiltonian for the 
KLM of the same form as (7). The present result is more 
general. For specific choices of PAM parameters, (7) re- 
duces to their result (and hence the eigenvalue spectra 
of two Hamiltonians become identical). The main point 
is that the equivalence in forms establishes a formal con- 
nection between the strong hybridization (mixed valence) 
and strong coupling regime of the PAM, and the strong 
coupling limit of the KLM. Given this connection, we rec- 
ognize that the bonding singlet, represented by a hole in 
H e ff, is a Kondo singlet in a mixed valence regime. 
It is interesting to remark that the large J/t limit of the 
KLM has instead been associated to the small U limit 
(weak coupling) of the PAM [15]. 

Sigrist et al [14] also proved that to leading order in 
t'i (1 < i < 3), the ground state of H e ff is ferromag- 
netic for any concentration of c electrons. To this end, 
they noticed that t[ is the only relevant term since it can 
change the spin configuration. The particular sign of t\ 
allowed them to apply the Perron- Frobenius theorem and 
demonstrate the ferromagnetic character of the ground 
state. Following exactly the same procedure, we can 
prove that the spin degeneracy of the ID PAM for t = 
and |e/| < |V| is lifted in a perturbative sense toward 
a ferromagnetic state for any concentration | < n < f 
with a total spin per site 2(n — i). 

The microscopic mechanism for the stabilization of the 
FM ground state is similar to the one found by Nagaoka 
[16] in the Hubbard model for dimension larger than one. 
Even though H e f f is a one dimensional model, there are 
two different ways to move one hole to a next-nearest 
neighbor site by hopping over one c-electron: either by 
two applications of t or by one application of t[. Only 
when the background is FM will both processes give rise 
to the same final state. If the sign of the final state is 
the same for both processes, the effective matrix element 
is reinforced (constructive interference) and kinetic en- 
ergy of the hole is lowered. This is indeed what happens 
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in H e ff due to the positive sign of t\. Therefore, the 
coherent propagation of the Kondo singlet (hole) is re- 
sponsible for the polarization of the spins which are not 
quenched by the conduction electrons. Since this coher- 
ent propagation is only possible for dense systems where 
the /-orbitals form a lattice, the stabilization of the FM 
state is excluded from any analysis which only considers 
the dilute limit (few /-orbitals). 

The same perturbative analysis used to prove the FM 
character of the ground state can also be used to show 
[14,17] that the low energy spin excitations of H e ff are 
described by a FM Heisenberg model with an effective 
exchange interaction: J = 2 ^[^- sin 2 (irp) — sin(2wp)], 
where p is the density of c-electrons. This effective 
Heisenberg model operates in a truncated Hilbert space 
of H e ff that only contains spin degrees of freedom. The 
separation between charge and spin degrees of freedom 
occurs because to leading order in t\, the eigenstates are 
still factorisable into their orbital and spin components. 

In summary, we proved that the PAM has a ferro- 
magnetic ground state in the strong hybridization and 
strong coupling regime for concentrations between half 
and three quarters filling. This wide region of concentra- 
tions clearly indicates that the microscopic mechanism is 
not related to an RKKY interaction as is expected for a 
strong mixed valence regime. Instead, the stabilization of 
the FM ground state is due to the coherent propagation 
of the Kondo singlets in a FM background. The mecha- 
nism is similar to the one operating in the Nagaoka [16] 
solution of the infinite U Hubbard model. 

It is important to ask what can we expect in higher di- 
mensions. In this case, the spin degeneracy of H(t = 0) 
is lifted to first order in t because there is no separation 
between the charge and the spin degrees of freedom. It 
is well known that the ferromagnetic Nagaoka [16] state 
is stabilized when one electron is added to the half filled 
system. Different works [18-20] indicate that a partially 
polarized FM state is favored for a finite range of con- 
centrations away from half filling. In addition, we have 
recently found numerical evidence of itinerant ferromag- 
netism in the PAM for the mixed valence regime and 
\V\<\t\ [5,6]. 

The current situation clearly merits reasking what are 
the possible scenarios when a dense system is away from 
half-filling to evolve from the localized to the mixed va- 
lence regime. In the localized regime, the hybridization 
can be considered as a perturbation and the correspond- 
ing perturbation theory tells us that the magnetic proper- 
ties are dominated by the RKKY interaction. When the 
system approaches the mixed valence regime, there are at 
least two different possibilities: a) A paramagnetic Fermi 
liquid, where the absence of magnetic ordering is just a 
consequence of the Pauli exclusion principle (only a small 
fraction of the /-moment is screened by the conduction 
electrons) [7] ; and b) A partially polarized FM metal (see 
also [5,6]) The existence of this second scenario as a pos- 
sible solution of the PAM provides a natural explanation 



for a number of U (US, USe, UTe [21], URu 2 _ a; M a; Si2 
with M = Re, Tc and Mn [22]) and Ce based (CeSi* [8], 
CeGe 2 [9], Ce(Rh 1 _ :c Ru ;c ) 3 B 2 [10], CeSii. 76 Cu . 2 4 [11], 
Ce 3 Bi4, and CeNio.sPto.2 [12]) materials which exhibit 
the coexistence of ferromagnetism and Kondo behavior. 

In addition, we have established a connection between 
the strong hyridization limit of the PAM and the strong 
coupling limit of the KLM which is valid in any dimen- 
sion. This connection provides a physical meaning for 
the strong coupling regime of the KLM. 
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